#ifndef __FDMATRIX_H__
#define __FDMATRIX_H__

void createFDmatrix(int iter_z,double dt,double dz, double material[][3],
               double phi[],double k0[],double Hpicar[],double H0[],double rhomx[][3], double A[], double b[],int lda);

#endif
